Theoretical investigations of a highly mismatched interface: the case of SiC/Si(001) 
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Using first principles, classical potentials and elasticity theory, we investigated the structure of a 
semiconductor/semiconductor interface with a high lattice mismatch, SiC/Si(001). Among several 
tested possible configurations, a heterostructure with (i) a misfit dislocation network pinned at 
the interface and (ii) reconstructed dislocation cores with a carbon sub-stoichiometry is found to 
be the most stable one. The importance of the slab approximation in first principles calculations 
is discussed and estimated by combining classical potential techniques and elasticity theory. For 
the most stable configuration, an estimate of the interface energy is given. Finally, the electronic 
structure is investigated and discussed in relation with the dislocation array structure. Interface 
states, localized in the heterostructure gap and located on dislocation cores, are identified. 

PACS numbers: 68.35-p, 81.05.Hd, 81.15.Aa 
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I. INTRODUCTION 



The misfit strain present in lattice mismatched epi- 
taxial layers has been widely studied because of its om- 
nipresence in the heterostructures used in device technol- 
ogy. Two goals are usually pursued, either to avoid strain 
for preparing long-lived devices, or to exploit the modi- 
fied electronic properties of strained layers for obtaining 
specific devices such as lasers.^' - The general picture of 
heteroepitaxy is well known. First, strained layers are 
grown on the substrate. Over a critical thickness, which 
depends on the elastic properties and the lattice mis- 
match between the two materials, plastic relaxation of 
the strain with formation of misfit dislocations becomes 
energetically favorable. In the case of small mismatch, 
the critical thickness can be large, with a low dislocation 
density. For example, a critical thickness of about lO'^ A 
is measured for a lattice mismatch of 1% in Ge^rSii-^ 
layers {x = 0.24)S, and the average separation between 
misfit edge dislocations is estimated to be about 390 A 
(for [001] layers). Each dislocation is far from the others 
and from the interface, and the system can be described 
by considering an ideal coherent interface, with strained 
layers. Hoekstra and Kohyama used this frame for in- 
vestigating the /3-SiC(001)/Al interface!^ In a scheme re- 
cently proposed, Benedek et al. introduced an additional 
correction for the effect of the misfit dislocations^, by 
comparing coherent and semi-coherent interfaces. How- 
ever, in the case of heteroepitaxy for largely mismatched 
systems, the critical thickness is very low and a dense 
network of misfit dislocations is present in the grown lay- 
ers. The interactions between dislocations and between 
the interface and the dislocations must then be explicitly 



taken into account. 

Despite the large lattice mismatch of ~20%, /3-SiC 
can be grown on Si(OOl), using different techniquesi^iii^i^ 
A recent High-Resolution Electron Microscopy (HREM) 
study of this system shows a locally abrupt inter- 
face, with the presence of a periodic array of misfit 
dislocations i^iSii^ These ones seem to be located directly 
at the interface, which points at a vanishing critical thick- 
ness for this heterostructure. It is difficult to gain addi- 
tional information from these experiments. In particu- 
lar, the atomic structure and the chemical environment 
at the interface, which deeply affect its physical proper- 
ties, still remain hardly accessible. Atomistic computa- 
tions are a possible solution for complementing the ex- 
periments. In addition, first principle calculations would 
allow the determination of the electronic properties of the 
heterostructure. In this respect, the SiC/Si interface may 
be considered as a model for a high lattice mismatch in- 
terface between covalent materials, just as silicon is usu- 
ally considered as the semiconductor prototype. Indeed, 
due to the peculiar lattice mismatch, almost equal to 
20%, the network dislocation pattern is extremely dense 
and can be modelled within a cell small enough to make 
an ab initio computation affordable. 

As far as we know, few studies have been devoted to 
the SiC/Si(001) system. Chirita et ali^ have investi- 
gated the strain relaxation and the thermal stability of 
the interface using the semi-empirical classical potential 
from Tersoffjii They have obtained a possible configura- 
tion of the atomic structure, stable up to 1000 K. How- 
ever, they assumed a stoichiometric interface, which may 
be a metastable state. Another study from Kitabatake 
with the same potential choice was exclusively focussed 
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on the first step of the formation of the interface, i.e. a 
pre-carbonization of the Si(OOl) surface prior to growthi^ 
A detailed investigation of the atomic structure and stoi- 
chiometry of the SiC/Si(001) interface with first principle 
techniques is still lacking. The large computational cell 
required for the modelling, notwithstanding the advanta- 
geous lattice mismatch, may explain the absence of such 
calculations. 

Recently, we have presented in a letter the most strik- 
ing results of the first ab initio study of the SiC/Si(001) 
intcrfacCi^^ Here, we provide a thorough description of 
our investigations, combining elasticity theory, classical 
potentials and first principles calculations. In particu- 
lar, after a complete explanation of the computational 
procedure, the full set of tested atomic configurations is 
presented. We have investigated the stoichiometry of the 
interface as well as its stability with few pseudomorphic 
SiC layers. The determination of the interface energy as 
well as the analysis of finite size effects in the slab ap- 
proximation are also presented. Finally, the computed 
electronic structure of the most stable configuration is 
discussed. 



II. COMPUTATIONAL PROCEDURE 
A. Model 

Considering that the lattice parameters of Si and 3C- 
SiC are respectively a(Si) = 5.4309 A and a(SiC) = 
4.3596 A3, the calculated misfit (a(Si) - a(SiC)) /a(Si) 
at the interface is 19.73%. Within the hypothesis of a co- 
herent interface, such a misfit is equivalent to a huge ten- 
sile stress in the SiC layers. A simple estimate shows that 
the elastic energy stored in the coherent film would be 
approximately 1.6 J/m^/layer.~ However, it is observed 
that the interface is semi-coherent, the stress being re- 
laxed by the introduction of a network of undissociated 
edge dislocations of Burgers vector b — a(SiC)/2(110), 
with dislocation lines lying along the [110] and [110] 
directions. ^'^'^ Due to the peculiar value of the misfit, 
approximately equal to ^ , the semi-coherence is obtained 
when the spacing between misfit dislocations is 56. The 
SiC/Si(001) interface is then modelled by matching a 
p(5 X 5)-SiC(001) slab (N layers) with a p(4x4)-Si(001) 
slab (M layers) along the [001] direction (z-axis), at the 
Si lattice parameter, optimized for bulk calculations with 
the chosen method. In the following, size conditions for a 
system will be simply designed by N/M. Periodic bound- 
ary conditions are applied along [110] and [110] directions 
(x- and y-axis) in any case. Atoms belonging to the top- 
most SiC and bottommost Si layers were fixed in bulk 
positions. However, in order to better relax the interface, 
the SiC topmost layer was allowed to move as a whole 
in all directions. This procedure was achieved at each 
step by (i) compute the average of atomic forces in this 
layer, (ii) apply this force to all layer atoms. Once the 
structures were converged, we also performed additional 



calculations to check the effect of the surfaces on the in- 
terface by releasing the constraints on both surfaces. 



B. Methods 

Classical molecular dynamics calculations were done 
with the semi-empirical Tersoff potentials'^, at the lat- 
tice parameter a(Si) = 5.43 A. This choice is currently 
the best available for describing solid SiC. A recent study 
has revealed that it reproduces DFT results on the for- 
mation energies and properties of native defects in SiC 
with a good accuracyJ^ Here, several system sizes were 
considered, because of the low CPU time and memory re- 
quirements for these calculations. The exploration of the 
configurational space was typically done with 12/8 and 
36/36 slabs, to ensure that the interaction between the 
interface and the surfaces is negligible. To estimate this 
interaction, sizes as large as 50/50 have been considered. 

First principles calculations were performecU^ within 
the Density Functional Theory (DFT) , in the Local Den- 
sity Approximation (LDA).'^ An energy cutoff of 40 
(160) Ry was used for the plane waves expansion of the 
wavefunctions (charge density). The reciprocal space in- 
tegration in the supercell Brillouin Zone (BZ) was done 
by considering only the F point. The core electrons were 
removed by using pseudopotentials, with s and p non- 
locality for Si and s non- locality for C— We used a cal- 
culated Si lattice parameter of 5.401 A, consistent with 
our ab initio pseudopotential. Surface atoms at both 
sides were saturated with hydrogens. The C-H and Si-H 
distances were optimized independently by fully relaxing 
surface saturated symmetric slabs of 11 layers, 8 atoms 
each. 5/5 and 7/7 system sizes were used, which corre- 
sponds to large scale calculations with at most 369 atoms. 
At variance with classical simulations, periodic boundary 
conditions were also applied along the [001] direction, a 
large vacuum space (9.1 A for 5/5, 8.2 A for 7/7) be- 
ing introduced to prevent spurious surface-surface inter- 
actions. All structures were considered converged when 
forces acting on atoms were less than 10"'' a.u. (0.005 
eV/A) and energy varied by less than 10~^ eV/atom. 



C. Configuration energy 

In this work, we have always compared the energies 
of systems with the same N/M. The numbers of C and 
Si atoms in the SiC slab may however change for two 
N/M slabs with different interface configurations. In ad- 
dition, the surface of the SiC slab may be either Si- or 
C-terminated. Hence, the energy Esiab of the N/M slab 
model described above, obtained from either methods, 
cannot be used directly for determining the most stable 
structure. Instead, a configuration energy Ea, given by 
the following expression, can be defined within the grand 
canonical frame: 
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III. STABILITY AND GEOMETRY 
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is the surface energy of the SiC part of the slab, Si- 
or C-terminated depending on the configuration at the 
interface and on the number of layers, whereas is the 
surface energy of the Si part, n^' and are the numbers 
of Si or C atoms in the SiC/Si slab. Finally, /x^' and ii'^ 
are the chemical potentials for each species. For the Si 
part of the slab, /i^' is equal to the bulk silicon chemical 
potential. Though and cannot be exactly deter- 
mined for the SiC part, it is possible to determine a range 
of thermodynamically allowed values;"^^ The following re- 
lations must be satisfied at the same time: 
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Here, the fi^ are the chemical potentials of the mono- 
elemental bulk phases, and are easily calculated. A7i/ 
is the heat of formation of silicon carbide: in this work, 
we used the experimental value, AHf = 0.72 eV. These 
two relations reveal that the Si and C chemical potential 
values are linked, and that the chemical potential of one 
species is always equal or lower than the chemical po- 
tential of the mono-elemental bulk phase. The allowed 
ranges are then 
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(/Xq' — AH/j^q) corresponds to C-rich (Si-poor) condi- 
tions and (^q', ^0 ~ ^^/) to C-poor (Si- rich) conditions. 
The values obtained in the present work for /Iq are listed 
in Table m here, for the ab initio results, we indicate the 
respective values as obtained from standard converged 
bulk calculations and as extrapolated from surface cal- 
culations with the scheme originally proposed by Fioren- 
tini and Methfessel.^° This will be useful to determine 
the different contributions in Eqn. (1) and the interface 
energy, as described below. It is worth noting that ex- 
treme caution must be adopted when estimating /Iq, as 
small (~ 1%) errors in this quantity are then multiplied 
by the usually large number of atoms at the interface, 
eventually producing non-physical results. The use of 
the linearized values of /ip in this work yields different 
configuration energy values with respect to our previ- 
ously published results-^^, without changing the original 
qualitative conclusions. 

Using the configuration energies allows to compare in- 
terface structures in the slab approximation, with differ- 
ent surface terminations and number of atoms. In the 
following, we discuss the stability of different configura- 
tions. 



In a previous work on the SiC/Si(001) interface, 
Chirita et al. proposed a possible geometry of the 
interface^", based on molecular dynamics calculations. In 
the choice of their initial atomic configuration, they have 
made several hypothesis. They assumed that the dislo- 
cations network is pinned directly at the interface, that 
the first (001) layer of the SiC slab is silicon-like, and fi- 
nally, that the interface remains perfectly stoichiometric. 
However, from HREM experiments, it is hard to extract 
such information. An interface with a pseudomorphic 
first layer followed by a misfit dislocations network is 
possible, as well as reconstructed misfit dislocation cores 
including a sub- or over-stoichiometry in carbon and/or 
silicon atoms. In our work, we have taken into account 
such possibilities, and our search for the most stable con- 
figuration has been broadened by considering a large set 
of initial geometries, by means of classical molecular dy- 
namics. First, only interfaces made with complete lay- 
ers were investigated. Fig. [T] shows 6 different possible 
starting geometries. The first ones, SI and S2, are stoi- 
chiometric configurations with the misfit dislocation ar- 
ray located directly at the interface. In SI, the first SiC 
layer is carbon-like whereas it is silicon-like in S2 (S2 is 
the initial configuration selected in the previous study^*'). 
In PI and P2, the misfit dislocation array is located in 
the second and third SiC layers, respectively. Those are 
obtained by removing one or two p(5x5) SiC layers, re- 
placed by p(4x4) pseudomorphic layers. Finally, we also 
compared the previous structures with configurations CI 
and C2, with a higher carbon concentration at the inter- 
face (carbonization). 

The configuration energy differences are reported in 
Table HIl Here, the calculated values are easily extracted 
from the classical dynamics, as the energy of each atom is 
a well defined quantity in the interaction potential model. 
The S2 configuration, investigated in a previous study 
is stable but is clearly not the lowest energy solution. 
Instead, the most stable structures are obtained for a 
stoichiometric system and a C layer at the interface, i.e. 
SI. For such structures, we found two different atomic 
configurations, represented in Fig. [2] It must be empha- 
sized that the two [110] and [lIO] directions, perpendicu- 
lar to the interface, are not equivalent due to the peculiar 
zincblende stacking. The first configuration, called Sla, 
has the lowest energy in the full range of chemical poten- 
tials. Along [110], bonds are formed between the extra 
C atoms located at the misfit dislocation core and Si 
atoms belonging to the second layer of the Si(OOl) slab. 
All these atoms are shifted toward the interface inducing 
large strains, especially in the softer silicon slab. Some of 
the Si atoms appear to be overcoordinated. Along [110], 
the main feature is a soft dimerization of the C atoms 
of the first SiC layer, on both sides of the misfit disloca- 
tion line. The C-C distance is shortened to 1.65-1.80 A 
which allows a reduction of the number of dangling bonds 
(DB). This interface reconstruction is nearly similar to 
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the model presented by Long et al.M- The second con- 
figuration, Sib, has a higher configuration energy, with 
6.55 eV in excess. The main difference with Sla concerns 
the [110] side. Here, no bonds are formed between the 
C atoms in the first SiC layer and the second layer of 
the silicon slab. Instead, the number of DBs belonging 
to C atoms is reduced by a stronger carbon dimeriza- 
tion along the [110] direction (clearly visible on the [110] 
view), where the C-C distance is reduced to 1.44-1.47 A. 
From Table |TT1 it seems that the occurrence of pseudo- 
morphic SiC layers before the introduction of misfit dis- 
locations is not favored, although the PI structure, with 
one pseudomorphic layer, is the third best solution. The 
segregation of carbon at the interface, tested with two 
configurations CI and C2, is associated with very high 
configuration energies. This is consistent with the gen- 
eral observation that the deposition of carbon on Si(OOl) 
surface does not lead to stable carbon layers but rather 
to the formation of silicon carbide by carbonization, the 
substrate supplying the Si atomsi^ 

Several information may be obtained from our results. 
First, one or several complete pseudomorphic SiC lay- 
ers are not favorable structures, which tends to indicate 
that the misfit dislocation network is located at the in- 
terface and no finite critical thickness can be defined. 
Secondly, all the favored configurations include a car- 
bon layer at the interface. This can be explained by 
the energy gained in the formation of SiC bonds, and 
the smaller carbon atomic radius with respect to silicon. 
Finally, we observed that the minimization of the config- 
uration energy is best realized by the reduction of DBs, 
especially those associated to carbon atoms. However, 
this DB saturation occurs via the formation of bonds be- 
tween initially remote atoms, and the shortening of their 
separation generates energetically expensive strain in the 
atomic structure. This is particularly true for the short- 
est and strongest carbon-carbon or silicon-carbon bonds. 

For the lowest energy configurations, Sla and Sib, the 
carbon atoms, located in the first SiC layer and in the 
vicinity of the misfit dislocation cores, present the most 
stretched bonds. We investigated the possibility of an 
energy lowering by removing these atoms. We initially 
considered the Sib configuration. A significant energy re- 
duction was obtained by removing successively C atoms 
along [110] (inside the ellipsis, in Fig. [21). The additional 
removing of the C-dimer row along [110] (the other el- 
lipsis) leads to an even more stable configuration, rep- 
resented in Fig. [31 The energy difference compared to 
Sla is -6.04 eV/ceU (-15.40 eV/cell) in C-rich (C-poor) 
conditions, respectively.^^ Starting from Sla, and remov- 
ing the C-row along [110] and the C-dimer row along 
[110] as described previously, the energy was also low- 
ered, and we obtained the same final structure. This 
configuration, which we called CSS for Carbon SubSto- 
ichiometric, shows remarkable features. Along the [110] 
direction, dangling bonds created by the removal of C 
atoms are eliminated with the formation of Si dimers 
2.48-2.54 A long. Compared to Sla, the Si atoms of the 



silicon slab, previously bonded to C atoms, recover bulk- 
like positions, thus minimizing the strain. Along [110], 
after the removal of the C-dimer row, the core of the mis- 
fit dislocation is made of 8- atoms rings, including seven 
Si atoms and a lone C atom. Almost all these atoms are 
fully coordinated, at the expenses of some bond stretch- 
ing, particularly on the Si atoms in the second layer of 
the SiC slab, and at the intersection of the two perpen- 
dicular dislocations. The Si-Si bonds range from 2.33 to 
2.62 A. We found that CSS is the most stable configura- 
tion using the Tersoff potential. Indeed, tests performed 
on structures with further C removal, or selected C/Si 
exchanges, brought no additional energy reduction. 

The Tersoff potential results were confirmed by the 
ab initio method. Owing to the huge computational ef- 
fort required to deal with the interface, our investigations 
were restricted to three configurations, i.e. the low energy 
stoichiometric Sla and Sib, and the most stable geom- 
etry CSS. First, we found that only Sla and CSS were 
stable, Sib relaxing spontaneously to Sla. It is known 
that classical potentials tend to stabilize a larger num- 
ber of metastable structures compared to first principles 
methodsj ^^'^'^ No major geometrical changes were found 
for both configurations relaxed with first principles, com- 
pared to the classical results. A better description of the 
atomic structure of the reconstructed dislocation cores 
was however obtained. Considering CSS and the [110] 
direction, the Si dimer lengths range between 2.42 and 
2.50 A. Along [110], the dislocation core is made of Si 
bonds 2.36-2.52 A long. Thus the ab initio relaxation 
yields a more compact 8-atom ring. The energy differ- 
ences between Sla and CSS are reported in Table IIIIl 
for both methods and different slab lengths. Despite the 
small sizes imposed by ab initio simulations, it appears 
that the energy difference is already well converged for a 
5/5 slab, with a 0.2 eV/cell variation between 5/5 and 
7/7. The ab initio calculation validates our primary re- 
sults, i.e. the CSS configuration is the most stable one, 
for the whole allowed range of the chemical potentials. 
Moreover, even in C-rich conditions, the energy differ- 
ence between the two geometries is much larger than the 
error associated with such computations. 

The analysis of the CSS configuration topology gives 
some hints for understanding this result. The first set of 
calculations has indeed shown that it is preferable to have 
C atoms at the interface in the first SiC layer. From the 
second set, it appears that misfit dislocation cores with 
only silicon atoms are favored, the stretching of carbon 
bonds being energetically expensive. In CSS, the first 
SiC layer is carbon-like, but there is only one C atom 
involved in the reconstructed dislocation cores. It also 
presents the peculiar characteristics that almost no atoms 
are sub- or over-coordinated, owing to the formation of a 
topological ring along one direction and the formation of 
a silicon dimer along the other. The CSS configuration is 
then the best candidate to represent the atomic structure 
of the SiC/Si(001) interface. 

A qualitative indication of the residual strain distribu- 
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tion at the interface can be obtained by inspection of the 
deformation of the converged system coordinates, with 
respect to the ideal bulk-Hke positions. In Fig. [H we 
report the layer puckering defined as the maximum de- 
viation in the direction perpendicular to the interface, 
as obtained from ab initio calculations for a 7/7 sys- 
tem: the major deformations are localized in Si, which 
has smaller elastic constants than SiC. The warping de- 
creases when moving aside the interface, in agreement 
with experimental resultsS^- on the structural characteri- 
zation of SiC films grown on a Si(OOl) substrate, which 
evidenced an internal roughness of individual SiC planes 
that diminishes away from the silicon substrate. Fig. [5] 
represents a comparison of the residual strain field at the 
interface for both the CSS and Sla geometries, evaluated 
in terms of atomic displacements from ideal bulk-like po- 
sitions. The superior efficiency of CSS in relieving the 
strain is clearly evidenced by simple inspection. 

Regarding the critical thickness, it is interesting to 
compare the CSS and PI configurations. Indeed, the SiC 
film in the CSS geometry is constituted by a substoichio- 
metric C layer, as if obtained from PI by removing four 
C atoms in a row. It is thus difficult to define precisely 
a finite critical thickness for this interface, characterized 
by largely reconstructed dislocation cores, although our 
results support the experimental findings of a dislocation 
network pinned at the interface!^ 



IV. EFFECT OF SLAB SIZE 

The large lattice misfit between /3-SiC and Si allows 
the investigation of the (001) interface with ab initio 
methods, the spacing between dislocations being about 
a(SiC)/2(110) ~ 15.4 A. There are only 4 (Si portion) or 
5 (SiC part) atoms per (110) edge. However, because the 
problem is two-dimensional, a slab layer includes 16 (Si 
part) or 25 (SiC part) atoms. Computationally, the num- 
ber of layers one can use to model the interface is then 
severely limited. Here, we managed to calculate at most 
a 7/7 interface, i.e. ~ 300 atoms. Such size is enough 
when considering a coherent interface. Here, the pres- 
ence of a periodic network of misfit edge dislocations at 
the interface induces a strain field in both the SiC and Si 
portions of the slab. It is usually assumed that the strain 
field penetrates each part by a distance of the order of 
the dislocation spacingi^ In our case, we should then in- 
vestigate a 14/11 slab, which is beyond the capabilities of 
available supercomputers. Our largest calculations per- 
formed on a 7/7 slab indeed reveal that the strain field 
is still not negligible at both ends of the slab: when al- 
lowing for a global relaxation, the two surfaces became 
slightly bent, due to residual strain, which may further- 
more be different for different core structures. The flat 
surface constrain adopted in our simulations is a valuable 
approximation to perform energy comparison: this con- 
strain however modifies the strain field of the dislocation 
arrays and a surface-interface interaction is present in the 



system. This interaction depends on the core structure 
at the interface, and may affect the relative stability of 
the configurations. 

In this part, we investigate the slab size effect by means 
of classical potential calculations and elasticity theory. 
The configuration energy, defined previously in Eqn. (1), 
may be written as the sum of four contributions: 

(h) = Ea + E, + Eel (h) + E,s (h) (6) 

Here, h is the slab size. Ea is a constant adhesive en- 
ergy between Si and SiC parts; Ec is the core energy 
of the misfit dislocation network; E^i is an elastic en- 
ergy due to the strain field; Eis is the interaction energy 
between surface and interface, which is zero for a slab 
including a large (infinite) number of layers. Only three 
terms depend on the slab size h. E^i is determined by 
using isotropic elasticity theory and a model of misfit 
dislocation arrays at the interface between a thin film of 
height h and an infinite substrate»2^ The dislocation core 
radius^^ is assumed to be equal to the Burgers vector 6, 
i.e. the in-plane SiC lattice parameter a(SiC)/2(110), in 
the present case. Provided that Ef-i is known for each 
h, Eis may be obtained from Eqn. [6] by calculating the 
configuration energy Ea with increasing slab size, for a 
chosen dislocation network. Here, the surface-interface 
interaction energy is determined separately for both the 
Si and SiC parts of the slab. We performed Tersoff po- 
tential calculations for slabs 31/n (increasing the Si part) 
and n/31 (increasing the SiC part), with n ranging from 
5 to 31, and for both the Sla and CSS configurations. 

The calculated interaction energies between surface 
and interface, Eis, are shown in Fig. [51 for both config- 
urations, in both portions of the slab, together with the 
sum of these terms. The SiC contribution appears larger 
than the Si one, for a given number of layers. Note that 
the interlayer spacing for SiC is about 20% lower than 
for Si, and more SiC layers are needed to get an equiva- 
lent contribution from the SiC and Si parts of the slab. 
For example, a 10/8 slab will have SiC and Si parts of 
nearly identical weights on Eis- A larger SiC contribu- 
tion is then expected for a lower number of layers. Here, 
the difference is important, in particular for CSS with an 
interaction energy almost three times larger for SiC than 
for Si. This could be explained by the geometry of CSS. 
The dislocation core along one direction is reconstructed 
with a ring of atoms, almost entirely located in the SiC 
part of the slab. The surface/interface interaction is then 
stronger in the SiC part. 

Our results confirm that the strain field penetrates by 
a distance of the order of the dislocation spacing^, and 
that ideally a 14/11 slab should be used. In fact, the sum 
of the SiC contribution for 14 layers with the Si contri- 
bution for 11 layers is less than 0.02 eV/cell, for both 
solutions. Assuming that the surface/interface interac- 
tion energy is mostly elastic, it is reasonable to consider 
that this quantity could be fairly estimated using classical 
potentials. Considering the total interaction energy Eis 
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in Fig. [6l for the Sla configuration, Eis is 0.66 eV/cell 
(2.36 eV/cell) for a 7/7 (5/5) slab, while for CSS, Ei, 
is 0.80 eV/cell (3.69 eV/cell) for a 7/7 (5/5) slab. In 
all cases, Eis is lower than the calculated configuration 
energy differences, and our previous conclusions on the 
stability of the CSS configuration remain valid. More- 
over, since the slab size effect is stronger for CSS than 
for Sla, using larger slabs will further increase the sta- 
bility of the CSS configuration. 

V. INTERFACE ENERGY 

As far as we know, there is no measured value of 
the SiC/Si(001) interface energy. Experimentally, it is 
possible to determine the bonding energy, which is re- 
lated to the interface energy, from wafer bonding exper- 
iments. However, a large range of values may be ob- 
tained, depending on the kind of SiC polytypes or sur- 
face terminations.— From the theoretical point of view, 
no value is available. Possible explanations are either 
the large size of the system that have to be dealt with 
first principles methods, or the difficulty to extract such 
energy from a slab calculation. A simplified frame for ob- 
taining the interface energy of mismatched interfaces has 
been proposed recentlyi^ The method however does not 
take into account a possible reconstruction or understo- 
ichiometry of the core of the misfit dislocations, which 
is mandatory for the SiC/Si(001) system, as discussed 
above. 

In order to extract the interface energy of the CSS 
configuration, the energies of the surfaces on both sides 
of the slabs have to be known. However, with ab ini- 
tio methods, only the total energy is accessible, and 
surface energies cannot be obtained directly. Instead, 
they were determined by extrapolating from slabs with 
an increasing number of layers, following the scheme 
proposed by Fiorentini and Methfessel;^ Si(001)-(lxl) 
and C:SiC(001)-(l X 1) hydrogenated surfaces were inves- 
tigated, since only the carbon terminated surface was 
relevant for the selected configurations. In order to min- 
imize computational errors, a large (c(4x4)) slab, with 
increasing thickness, from 7 up to 15 layers, was used: 
the extrapolated values for the chemical potentials are 
close to those obtained from bulk calculations for Si and 
SiC respectively (see Table IJ). To be consistent with 
the SiC/Si interface slab calculation, each surface atom 
was saturated by 2 symmetric hydrogens. As a conse- 
quence, the surface energies reported in Table HVl include 
the energies of the pseudo- hydrogens. We found a C-H 
distance of 1.09 A and a bond angle HCH of 100.6° for 
the C:SiC(001)-(lxl), and a Si-H distance of 1.47 A and 
a bond angle HSiH of 101.4° for the Si(001)-(1 x 1). We 
point out that we applied constraints in order to keep a 
symmetric dihydride Si(001)-(lxl) surface, although the 
most stable configuration is cantedi^ We forced the sym- 
metric geometry to quickly recover a bulk-like behavior, 
as required by the small slab sizes adopted in the cal- 



culations, since the correct canted configuration extends 
deeply in inner layers. 

The interface energy ai is obtained as the converging 
value of the configuration energy for large N/N slabs. In 
the classical case, it is easily estimated with large slabs, 
where the interaction between surface and interface is 
negligible, and all the elastic energy can be considered 
fully pertinent to the inner layers of the slab. In this 
case, we obtain ct/ ~ 22.7 eV/cell, (i.e. aj ~ 1.6 J/m^). 

Using ab initio methods, we are limited to small slabs, 
5/5 and 7/7. In principle, the determination could be 
done in a similar way than for surface energies, using an 
extrapolating techniquci^ However, here, the available 
slab sizes are too small for that purpose. In fact, both 
elastic and interface-surface energies considerably change 
between the 5/5 and the 7/7 slabs. As a consequence, a 
linear regression would yield misleading values for the 
slope and intercept constants. The extrapolating scheme 
could be used for slabs large enough to have constant 
elastic and surface-interface interaction energies, i.e. for 
N/N slabs with N greater than 14 for example. Instead, 
we directly calculate the interface energy for a 7/7 slab, 
using the chemical potentials obtained via linear extrapo- 
lation, as discussed above. We computed aj ~ 23.0 (22.6) 
eV/cell in C-rich (C-poor) conditions. This quantity is 
not exactly the interface energy. Indeed, for a 7/7 slab, 
the surface-interface interaction is not negligible and the 
strain energy is still not fully converged. These contribu- 
tions have been determined in the previous section: for 
a 7/7 slab in CSS geometry the excess surface- interface 
interaction is 0.8 eV/cell (see Fig. 4), and the calculated 
residual elastic energy between an infinite interface and 
a 7/7 slab is 0.2 eV/cell. The corrected interface energy 
is then ai ~ 22.4 (22.0) eV/cell (i.e. 1.58 (1.55) J/m^) in 
C-rich (C-poor) conditions. The agreement between clas- 
sical and ab initio values is surprinsigly good and maybe 
fortuitous given the technical difficulties associated with 
the interface energy determination in the ab initio case, 
or the use of the Tersoff potential that does not include 
properly the electronic contributions. 

We wish to stress again at this point that the ab initio 
values provided for the interface energy are mere esti- 
mates, as small variations in e.g. the chemical potentials 
can induce large errors and even non-physical negative 
interface energies. Our choice in presenting these results 
has been to keep consistency between the chemical poten- 
tials used to eliminate the surface and bulk contributions: 
for this reason, we always adopted the values for /ig as 
obtained in the linear extrapolation scheme^^, for large 
supercells, that allowed a good fcy sampling. 



VI. ELECTRONIC PROPERTIES 

We now move to the description of the electronic struc- 
ture of the mismatched heterostructure. The interface 
configuration determines the electronic properties of the 
system: the presence of defects such as misfit disloca- 
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tions can induce interface states in the band gap, that 
may severely modify device performances. Indeed, our 
results indicate that a number of interface derived states 
lay in the forbidden energy gap for the most stable dis- 
location network, although the number of DBs is in this 
case minimized. 

For a 5/5 slab, the valence band widths (VBW) of Si- 
and SiC-derived bulk states compare fairly well for frozen 
and free surfaces, although they underestimate the re- 
spective bulk calculations (Table |V|) . Increasing the slab 
thickness to 7/7 leads to VBW variations of only ~ 3%. 
We were able to estimate an error due to the slab approx- 
imation of ± 0.3 eV on LDA eigenvalues. A 7/7 slab is 
then large enough to get a good description of both the 
Si and SiC part, and the perturbation induced on the 
electronic structure by the interface configuration. In 
Fig. [71 the spatially projected Density Of States (DOS) 
at the interface is compared with the DOS obtained in 
inner layers at the Si and SiC sides of the slab. Several 
states lay in the band gap of the heterostructure, as ob- 
tained by alignment of the respective bulk valence bands 
(see Table |V] and the peak above the valence band top 
at r, highlighted by the arrow in Fig. [7]). The Highest 
Occupied (HO) and Lowest Unoccupied (LU) states are 
located at 0.7 and 1.1 eV above the valence band top at 
r. These two states are also localized in the core of the 
edge dislocations, as a result of the large difference in 
charge transfer between Si-C and Si-Si bonds. In Fig. [51 
the charge density plot of the HO state along the two dis- 
locations directions is represented. The charge density is 
mainly localized on atoms of the [110] core dislocation, 
while no density is observed around the other core in the 
perpendicular direction. The opposite situation is found 
in other bonding states localized in the system forbid- 
den gap and in particular for the LU state (see Fig. ^ 
whose charge density pertains to the dislocation core lay- 
ing along the [110] direction. Clearly, these states are 
true interface states, resulting from the reconstructed dis- 
location cores. It is worth noting that they would not be 
obtained from a coherent interface calculation, and that 
the determination of the electronic structure requires the 
atomic characterization of the misfit dislocation network. 

It has been recently pointed out that a peak is observed 
in EELS and XPS around 0.8 eV.^^ Our results on the 
presence of interface states in the heterostructure forbid- 
den gap may give an explanation for this experimental 
evidence. 



VII. CONCLUSION 

In conclusion, we have characterized the energetics, 
the atomic and electronic structure of the SiC/Si(001) 
interface. We performed first principles DFT-LDA, clas- 
sical potential and elasticity theory calculations. The 
most stable atomic configuration is in agreement with 
experimentsi'^'^, where an array of misfit dislocations 
pinned at the interface is observed. Actually, from our 



results, we also predict that the dislocation core is char- 
acterized by a substoichiometry in carbon atoms. Addi- 
tional experiments are needed to confirm our proposed 
configuration: our relaxed structures may be used as in- 
put for simulating HREM experimental images and com- 
plement the results. 

We furthermore estimated the interface energy from 
ab initio calculations, for a non-coherent interface, in the 
case of a multicomponent system, like SiC: this quantity, 
which is hardly accessible from experiment, has been here 
evaluated for the SiC/Si(001) interface, although with 
a non-negligible uncertainty. This is to our knowledge 
the first ab initio determination of interface energies at a 
mismatched semiconductor heterostructure. 

Several electronic interface states, calculated at F, have 
been identified. These states, located on the core of the 
misfit dislocations, may infiuence the electronic and op- 
tical properties of the interface. 

We point out that the approach we used in this study, 
i.e. the combination of elasticity theory, classical poten- 
tial, and ab initio methods, may be easily adapted for 
other systems of interest, in particular for systems with 
large mismatch, where a coherent interface model is not 
suited, or for semiconducting systems, where core recon- 
structions are expected. 
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Table captions 



TABLE L Optimized chemical potentials and lattice con- 
stants for Si and 3C-SiC bulks, computed with the Tersoff 
potential and the ab initio method. For ab initio, the bulk cal- 
culations (column "bulk") and the value obtained from linear 
extrapolation'^" (column "Jin") are indicated for comparison. 



TABLE II: Configuration energy differences, from classical 
dynamics, and variation of C and Si atoms, for the first set 
of configurations and a 12/8 system (see Fig. [!}, with the 
geometry Sla taken as the reference. 



TABLE III: Configuration energy differences for CSS com- 
pared to Sla, within different conditions, and slab size. As 
indicated in the text, numbers for the classical simulations 
are slightly different from those reported in Ref.— , because 
of improved convergence. The ab initio values are obtained 
with the linear extrapolation method^S^ and the consistently 
derived chemical potentials were used. 



TABLE IV: Surface energies for the systems used in the 
slab calculations, computed with the ab initio method. The 
pseudo-hydrogen energy is included, as detailed in the text. 



TABLE V: Valence band Width (VBW) and energy gap 
(Agap) of the relaxed structures in eV. 





Torsoff 


Ab initio 






fio (eV) ao (A) 


/io(bulk) (Ha) ^io(lin) (Ha) 


ao (A) 


Si 


-4.630 5.432 


-3.96611 -3.96877 


5.401 


SiC 


-12.374 4.318 


-9.64469 -9.65374 


4.334 



Table H 



Config. 






AEa (eV/cell) 
C-rich C-poor 


Sla 








0.0 0.0 


Sib 








6.55 6.55 


PI 





-9 


16.40 9.92 


P2 


-9 


-9 


40.09 40.09 


S2 








37.91 37.91 


C2 


-25 


16 


63.74 93.26 


CI 


-25 


25 


63.88 99.88 



Table HI 





Tersoff 


Ab initio 


(eV/cell) 


12/8 36/36 


5/5 7/7 


C-rich 


-6.04 -6.03 


-5.90 -5.69 


C-poor 


-15.40 -15.39 


-15.26 -15.05 



Table HIl 



Es (Ha/at-Ha) 


C-rich C-poor 


Si(001)-(lxl) 
C:SiC(001)-(lxl) 


-1.12208 
-1.14212 -1.12890 



Table |TV] 



System 


VBW (eV) 
SiC Si Total 


A9ap(eV) 


Si bulk 


- 12.08 - 


- 


SiC bulk 


15.68 - 




5/5 frozen surf. 


14.67 10.91 15.31 


0.50 


5/5 free surf. 


14.65 9.95 14.68 


0.96 


7/7 frozen surf. 


15.06 11.25 15.57 


0.40 



Table |V] 
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Figure captions 



FIG. 1: First set of initial configurations. White (black) cir- 
cles represent the Si (C) atoms. The dashed line shows the 
location of the interface. 



FIG. 8: Isosurface (medium grey) of the highest occupied 
state projected along [110] (left) and [110] (right). Black 
(light grey) spheres indicate C (Si) species. Surface atoms 
are saturated with hydrogens (small white spheres). 



FIG. 2: [110] (left) and [110] (right) side views of three re- 
laxed SiC/Si(001) interface configurations. Light grey (black) 
spheres show silieou (carbon) atoms. The dashed thick lines 
mark the location of the extra atomic planes introduced by 
the misfit edge dislocations. Note that the represented bonds 
are drawn solely on the basis of a distance criterion and are 
not necessarily indicative of a true chemical bond. Removal of 
the C atoms inside the ellipses leads to the CSS configuration. 



FIG. 3: [110] (top) and [110] (bottom) perspective side views 
of the most stable CSS SiC/Si(001) interface configuration. 
Light grey (black) spheres show silicon (carbon) atoms. 



FIG. 4: Layer warping for a 7/7 SiC/Si(001) interface. In 
abscissa, the mean value of the coordinate normal to the in- 
terface plane for a given atomic layer is indicated. 



FIG. 5: Contour plot of the strain field at the 7/7 SiC/Si(001) 
interface projected along the [110] direction for both the CSS 
(top) and Sla (bottom) geometries, evaluated in terms of 
atomic displacements from ideal bulk-like positions. To en- 
hance comparison, the displacement integrated over a super- 
cell plane of fixed height from the interface is indicated. Two 
lateral replicas are indicated for both configurations, with the 
SiC part higher and Si lower; brighter regions depict larger 
distortions: they occur slightly below the interface, at the Si 
first layer. 



FIG. 9: Isosurface (medium grey) of the lowest unoccupied 
state projected along [110] (left) and [110] (right). Black (light 
grey) spheres indicate C (Si) species. Surface atoms are sat- 
urated with hydrogens (small white spheres). 



FIG. 6: Interaction energy between surface and interface as 
a function of the slab size, for the Sla (left) and CSS (right) 
configurations. 



FIG. 7: Calculated DOS for the CSS configuration and a 7/7 
slab, projected on the Si and SiC layers at the interface (top 
panel), compared with the DOS for bulk SiC and Si (middle 
and bottom panels). The width of Si and SiC bands, as well 
as the energies of HO and LU states, are shown in the figure. 
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Figure [5] 
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